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The transverse momentum and pseudorapidity dependences of partonic ' Mach-like' shock wave 
have been studied by using a multi-phase transport model (AMPT) with both partonic and hadronic 
interactions. The splitting parameter D, i.e. half distance between two splitting peaks on away side 
in di-hadron azimuthal angular (A</>) correlations, slightly increases with the transverse momentum 
of associated hadrons (py ssoc ), which is consistent with preliminary experimental trend, owing to 
different interaction- lengths/numbers between wave partons and medium in strong parton cascade. 
On the other hand, the splitting parameter D as a function of pseudorapidity of associated hadrons 
(ij° ss ™), keeps fiat in mid- pseudorapidity region and rapidly drops in high-pseudorapidity region, 
which is as a result of different violent degrees of jet-medium interactions in the medium that has 
different energy densities in the longitudinal direction. It is proposed that the research on the 
properties of ' Mach-like' correlation can shed light on the knowledge of both partonic and hadronic 
interactions at RHIC. 

PACS numbers: 12.38.Mh, ll.10.Wx, 25.75.Dw 



An exotic matter predicted by Quantum Chrome-Dy- 
namics (QCD) [l[ may be created in the early stage of 
central Au+Au collisions at ^smn = 200 GeV at the Rel- 
ativistic Heavy-Ion Collider (RHIC) at the Brookhaven 
National Laboratory 0j. The sufficient experimental ev- 
idences 

HHHH support that the new matter is not an 
ideal gas but a strongly-interacting dense partonic mat- 
ter (named as sQGP) under extreme temperature and 
energy density. 

Jet has been proved to be a good probe to investigate 
the characters of new matter in RHIC experiments @. 
Calculations based on pQCD predicted that high energy 
partons (jets) traversing a dense QCD medium lose en- 
ergy through induced gluon radiation 0, where the en- 
ergy loss is expected to depend strongly on the color 
charge density of the created system and the traversed 
path length of the propagating jets. The energy loss phe- 
nomenon (i.e. jet quenching) has been observed by some 
experimental probes. In di-jet azimuthal angular (A</>) 
correlations, one back-to-back hard jet disappears due to 
large energy loss in the dense medium Q . According to 
the conservation law of energy, the lost energy should 
be redistributed in the medium 0, [n| HH, El, which 
has been reconstructed by di-hadron azimuthal angular 
correlations of charged particles in experiments [l3| . Re- 
cently, an interesting ' Mach-like' structure (the splitting 
of the away side peak in di-hadron A(f> correlation) has 
been observed in both di-hadron and three-hadron az- 
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imuthal angular correlations in central Au + Au colli- 
sions at ^jvw = 200 GeV 0GSG1, 

which is attracting 

many theorists' attentions. 

Recently considerable theoretical efforts have been put 
into the topic 0, E, H US El- Among these re- 
searches, the transverse momentum (pr) and pseudora- 
pidity (77) dependences of 'Mach-like' structure are two 
subjects of much concern. Cherenkov-like gluon radi- 
ation model has been suggested to produce a conical 
structure [IU . It gave the dispersion relations for differ- 
ent choices of masses of two massive scalars coupled to a 
massless scalar. The dispersion relations are space-like in 
low momentum and approaches the light cone as momen- 
tum is increased. With the dispersion relations, the cone 
angle with respect to the jet direction becomes narrow 
rapidly with the increasing of the momentum of the emit- 
ted particle. The medium-induced gluon bremsstrahlung 
radiation H^j has been used to study the di-hadron 
correlations. In Ref the large angle medium- induced 
gluon radiation, gave a picture of the radiative angle rel- 
ative to the jet axis becoming narrow with the increasing 
of radiative gluon energy, but did not show a two-peak 
structure on away side. In Ref (20| , a simple generaliza- 
tion of Sudakov form factor to opaque media presented 
non-Gaussian shapes of away-side in di-hadron azimuthal 
correlations and a centrality dependence of the position 
of splitting peaks which can match experimental data 
well. However, there are no discussions about the px 
dependence of non-Gaussian shape in the model. In 
a shock-wave model proposed by Shuryak and his col- 
laborators linearized hydrodynamic equations were 
solved in a region far from the jet where the perturba- 
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tions and gradients are small. In general, the strength of 
the diffusion mode relative to the sound mode is directly 
proportional to the entropy produced by jet-medium in- 
teractions. Only when there is no significant entropy pro- 
duced by strong jet-medium interactions, hydrodynamic 
fields can be induced and peaks at the Mach angle in di- 
hadron azimuthal angular correlations can be revealed. 
It is noticeable that the results are very sensitive to the 
parameters chosen, especially need very large values for 
the energy loss to produce the correlation functions simi- 
lar to the experimental results. The position of the peak 
on away side, which corresponds to splitting parameter 
D (the half distance between two peaks on away side in 
a di-hadron azimuthal angular correlation), is related to 
the speed of sound c s . The value of c s used in the calcu- 
lation is a time-weighted average (c s ~ 0.33) consistent 
with the expected value for RHIC. As a result, the posi- 
tions of the peaks at away side are located at Acj> = n ± 
arccos(c s ) « 1.9, 4.3 (rad) and nopx dependence. On the 
other hand, a probable image of ' Mach-cone' shock wave 
was given as a function of rapidity [22j|, which implied 
that the splitting parameter D should become smaller 
with the increasing of the rapidity (y), arising from a 
folding of condition possibility p(y) and rapidity struc- 
ture of 1 Mach-cone' correlation. In our previous works, 
we have shown a partonic ' Mach-like' shock wave can be 
produced and evolved in strong parton cascade, devel- 
oped by hadronic rescattering, in two- and three-particle 
azimuthal correlation analysis [23|, [U Hf| ■ In this work, 
we report a study of px and r\ dependences of partonic 
1 'Mach-like' shock waves by using A Multi-Phase Trans- 
port model (AMPT) [2(|. The p T dependence of split- 
ting parameter (D) that D increases with p^? so slightly 
is attributed to different interaction-lengths/numbers be- 
tween wave partons and medium in strong parton cas- 
cade and r\ dependence of splitting parameter (D) that 
D keeps flat in mid-pscudorapidity and drops quickly in 
high-pscudorapidity is due to different strengths of jet- 
medium interactions in the medium which has different 
energy densities in the longitudinal direction. 

AMPT model [2(| is a hybrid model which consists 
of four main components: the initial conditions, par- 
tonic interactions, the conversion from partonic matter 
to hadronic matter and hadronic rescattering. The ini- 
tial conditions, which include the spatial and momen- 
tum distributions of minijet partons and soft string ex- 
citations, are obtained from the HIJING model [27l ]. 
Excitations of strings melt strings into partons in the 
AMPT version with string melting mechanism [28| (abbr. 
l the melting AMPT model'). Scatterings among par- 
tons are modelled by Zhang's Parton Cascade model 
(ZPC) [29], which at present only includes two-body scat- 
terings with cross sections obtained from the pQCD cal- 
culations with screening mass. In the default version 
of AMPT model [13|(abbr. 'the default AMPT model), 
minijet partons are recombined with their parent strings 
when they stop interactions, and the resulting strings are 
converted to hadrons by using the Lund string fragmen- 



tation model [3l[ . In the melting AMPT model, a simple 
quark coalescence model is used to combine all partons 
into hadrons. Dynamics of the subsequent hadronic mat- 
ter is then described by A Relativistic Transport (ART) 
model [l^l ■ Details of the AMPT model can be found in 
a recent review (26| . Previous studies have shown that 
the partonic effect could not be neglected and the melt- 
ing AMPT model is much more appropriate than the 
default AMPT model when the energy density is much 
higher than the critical density for the predicted phase 
transition [26|, [28|, [33j| . In the present work, the parton 
interaction cross section in the melting AMPT model is 
assumed to be 10 mb. 

Our analysis in this paper are based on di-hadron az- 
imuthal angular correlations between a high px hadron 
(trigger hadron) and low px ones (associated hadrons). 
The analysis method is similar to that the experimenters 
did [HI, [HI]. The pairs of associated and trigger par- 
ticles in the same events were accumulated to obtain 
Acf> = 4>assoc — 4>trig distributions. Mixing-event tech- 
nique was applied to construct background which is ex- 
pected mainly from elliptic flow [l3|, ■ In this method, 
we mixed two events which have very close centralitics 
into a new mixing event, and obtained A<f> distribu- 
tion which is regarded as the corresponding background. 
A Zero Yield At Minimum (ZYAM) assumption was 
adopted to subtract the background as did in experimen- 
tal analysis [l5j| (See Ref. (23| for our detailed analysis 
techniques). In this work, we took two sets of AMPT 
models (the default and melting AMPT models) to sim- 
ulate the central Au+Au collisions (0-10%) at ^/s^n = 
200 GcV. The px window cuts for trigger and associated 
particles are 3 < p% i9 < 4 GeV/c and < p T ssoc < 3 
GeV/c within pseudorapidity window |^ tri 9^ ssoc | < 1 for 
Px dependence analysis, and 2.5 < p^ lg < 6 GeV/c and 
1 < p9j? soc < 2.5 GeV/c within pseudorapidity window 
|t7*™ 9 | < 1 for rj assoc dependence analysis. 

Figure [1] shows p T ssoc dependences of splitting parame- 
ter D in Au + Au collisions (0-10%) at V^vw = 200 GeV 
in the melting and default AMPT model with/without 
hadronic rescattering (i.e. ART process). Here the trig- 
ger and associated particles are with 3 < p*^ 19 < 4 
GeV/c and < p T ssoc < 3 GeV/c respectively, both 
within pseudorapidity window |jy*"9>o ssoc | < 1. There 
is no significant difference of splitting parameter D be- 
tween the cases with and without hadronic rescattering 
in the melting AMPT model, but D in the default AMPT 
model is smaller than that in the melting AMPT model. 
It is consistent with our previous results (23| in which 
hadronic rescattering alone gives much smaller splitting 
amplitudes than that in the melting AMPT model which 
shows good D values matching experimental data. Here 
the melting AMPT model gives much better splitting pa- 
rameters D which are consistent with the preliminary ex- 
perimental trend (34| , but different from the calculations 
of Casalderrey-Solana et al. [TtJ and Koch et al. [Hj]. It 
is noticeable that the p T ssoc dependence of splitting pa- 
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FIG. 1: (Color online) The splitting parameter D as a func- 
tion of p^ aoc for trigger hadrons of 3 < p% iB < 4GeV/c 
from AMPT simulations for Au + Au collisions (0-10%) at 
y/sj^N = 200 GeV. Full circles: the melting AMPT model 
with hadronic rescattering (i.e. ART process); open circles: 
the melting AMPT model without hadronic rescattering; full 
squares: the default AMPT model with hadronic rescattering; 
band: preliminary experimental trend [34| . 
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2: (Color online) AMPT model calculations of < 
possoc ^ versus Acj> for away side in di-hadron A<f> corre- 
lations with pr window cut ( 3 < p% %g < 4 GeV/c and 
< p^ ssoc < 3 GeV/c) in Au + Au collisions (0-10%) at 
y/s^F = 200 GeV. Full circles: the melting AMPT model 
with hadronic rescattering ; open circles: the melting AMPT 
model without hadronic rescattering; full triangles: the melt- 
ing AMPT model without parton cascade (i.e. ZPC process) 
and hadronic rescattering; full squares: the default AMPT 
model with hadronic rescattering; band: preliminary experi- 
mental trend [34l |. 

rameter D is sensitive to parton cascade (Note that it 
was found that more long-lived partonic phase and big- 
ger parton interaction cross section can give more bigger 
D values in our previous work [H]].), but not to hadronic 
rescattering in the melting AMPT model. Therefore, the 
relation of D vs. p9j? soc can be a good probe that not 



only distinguishes different physical scenarios but also 
sheds light on the properties of partonic interactions in 
the early stage of the rclativistic heavy ion collisions. 

Figure [2] presents the relation of average transverse 
momentum of associated particles (< p^ ssoc > = 

^'^s" s ^° cPT ' ass ° c ) versus A0 for away side. There is a 
distinct dip structure at A</> = it for all model condi- 
tions except the melting AMPT model without parton 
cascade and hadronic rescattering, which indicates that 
away jets arc depleted by strong parton cascade in the 
melting AMPT model (or by hadronic rescattering in the 
default model) and its lost energy is expected to excite 
the formation of partonic (or hadronic) ' Mach-like' shock 
waves [25j . It is apparent that hadronic rescattering can 
soften the associated hadrons by comparing the results 
between with and without hadronic rescattering in the 
melting AMPT model. The depletion structure implies 
that harder associated hadrons prefer larger angles with 
respect to the away-jet direction, which results in the in- 
creasing trend of splitting parameter vs p9j? soc as shown in 
Figure [TJ In order to explain the dip structure, we make 
a simple MC simulation. Figure [3] gives the initial parton 
number density which is projected into x — y plane in the 
melting AMPT model. Once shock wave is generated at 
some points p(xo,yo), the wave- front direction would be 
symmetrical with respect to the away-jet direction (at 
Acf) = it rad). The average interaction-lcngths/numbcrs 
of wave-front partons interacting with the medium are 
with maximums at A<p = tt rad (i.e. in the away-jet di- 
rection), as shown in panel (a) and (b) of FigureHJ which 
is the cause of the dip structure of < p^? soc > vs A</> for 
away side. Though < p%? soc > vs. A</> in the final state 
is a result of both parton cascade and hadronic rescat- 
tering in the melting AMPT model, it is interesting that 
a combination research between < p?j? soc > vs A0 and D 
vs. p9j? soc may provide a potential tool to extract the re- 
spective information from two different interaction levels 
(i.e. partonic and hadronic stages). 

As we have already shown that the origin of partonic 
1 Mach-like' shock waves could be attributed to the big 
partonic interaction cross section which can couple par- 
tons together to show a hydro-like collective behavior in 
our previous papers (24l . t25| . and our used model is a 
MC dynamical partonic transport model which can har- 
moniously work in the regions close to and far from jet, 
whereas linear hydrodynamical approximation is only ap- 
plicable in the region far from jet as discussed above. 
Therefore, it is not strange that such an increasing de- 
pendence, which is different from the predicted by a hy- 
drodynamic calculation [17| . is presented in the present 
model. 

In the above studies on pt dependences of D, both the 
trigger and associated hadron both are located at mid- 



pseudorapidity (\rj 



trig,assoc I 



< 1). When one scans the 



whole pscudorapidity region to look for all hadrons which 
are associated with a mid-pseudorapidity trigger hadron, 
how about the associated hadrons in other pscudorapid- 
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FIG. 3: (Color online) The initial distribution of parton num- 
ber density in mid-rapidity for Au + Au collisions (0-10%) at 
^/s^Tzv = 200 GeV in the melting AMPT model, which has 
been projected into x-y plane, where p(a;o,t/o) is the genera- 
tion point of Mach-like shock wave and R is the scale size of 
the system. 
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FIG. 4: (a): The average interaction-lengths of a wave parton 
undergoing the medium in different A(f>; (b): The average 
interaction-numbers of a wave parton undergoing the medium 
in different A(j> for Au + Au collisions (0-10%) at y^ijvjv = 
200 GeV in the melting AMPT model. 
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FIG. 5: (Color online) The AMPT model calculations in 
Au + Au collisions (0-10%) at ^/S^" = 200 GeV. (model 
cuts: 2.5 < v% la < 6 GeV/c, 1 < p£ ssoc < 2.5 GeV/c, 
and \rj trlg \ < 1). (a): the pseudorapidity dependences of 
splitting parameter D; (b): the pseudorapidity distributions 
of associated hadron yields for away side. Full circles: the 
melting AMPT model with hadronic rescattering (i.e. ART 
process); open circles: the melting AMPT model without 
hadronic rescattering; full square: the default AMPT model 
with hadronic rescattering; lines: the corresponding Gaussian 
fitting functions. 

ity regions? Are they waked by an away-jet? Figurc[5](a) 
shows the dependences of splitting parameter D on the 
pseudorapidity of associated hadrons (7y assoc ) ; for a given 
trigger hadron within l^*™ 9 ! < 1. The splitting parame- 
ter D is almost unchanged in a wide mid-pseudorapidity 
region, but goes down quikly at high absolute value of 
pseudorapidity. For the splitting parameter D in the de- 
fault AMPT model with hadronic rescattering, it is lower 
than that in the melting AMPT model and not compa- 
rable with experimental results, which further demon- 
strates that strong parton cascade process is necessary 
for reproducing reasonable splitting parameters. Fig- 
ure [5] (b) gives pseudorapidity distributions of associated 
hadron yields for away side with a trigger hadron in mid- 
pscudorapidity. The associated yield on away side in the 
melting AMPT model is higher than that in the default 
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FIG. 6: AMPT model calculations in Au + Au collisions 
(0-10%) at ^/sjvjv = 200 GeV. (a): conditional probability 
P{ii) to find the away side associated parton (1 < p^, ssoc < 
2.5GeV/c) at pseudorapidity rj if the trigger parton (2.5 < 
p-p 9 < QGeV/c) is found at \rf ria \ < 1 in the initial state of 
melting AMPT model; (b) : the pseudorapidity distribution of 
initial partons (hadrons) (0 < pr < 2.5GeV/c) before parton 
cascade (hadronic rescattering). solid: 77 distribution of initial 
partons in the melting AMPT model; dash: 77 distribution of 
initial hadrons in the default AMPT model. 

AMPT model, but there is no distinct difference between 
the two cases in the melting AMPT model with/without 
hadronic rescattering. It is interesting that the associ- 
ated away-side yields in the melting AMPT model can 
be extended to high pseudorapidity region (\rj\ > 4), but 
those in the default AMPT model are confined within the 
region \rj\ < 4, which indicates that the medium at high 
pseudorapidity can be waked by an away-side jet due to 
strong parton interactions. 

In panel (a) of Figure O it shows the conditional 
probability P(rj) to find the away side associated par- 
ton (1 < p9f soc < 2.5GeV/c) at some pseudorapidity 
77, if a trigger parton (2.5 < < 6GeV/c) is found 



at |t7*" 9 | < 1 in the initial state before parton cas- 
cade, which is consistent with the LO pQCD calcula- 
tions _22j. It means that the away-jet is with a certain 
possibility to be in non-mid-pseudorapidity. At the same 
time, the pseudorapidity (77) distribution of initial par- 
tons (hadrons) (0 < pr < 2.5GeV/c) before parton cas- 
cade (hadronic rescattering) is also presented in panel 
(b) of Figure [SJ Since 77 distribution of particle number 
is proportional to the 77 distribution of energy density, 
the interactions between away side jets and associated 
particles should be more violent in mid-pseudorapi dity 
region than in high \rj\ region. In our previous work [25j, 
the formation of partonic ' Mach-like' shock waves stems 
from strong parton cascade, therefore the rj assoc depen- 
dence of splitting parameter is expected to result from 
different violent degrees of jet-medium interactions in the 
medium which has different energy densities in the longi- 
tudinal direction. From Figure 0(a) and (b), the rj assoc 
dependence of the splitting prameter can identify the par- 
tonic or hadronic origin of the Mach-like behavior, and 
provide the information about minijet production mech- 
anism and the formed matter density distribution in the 
longitudinal direction. 

In conclusion, the properties of partonic ' Mach-like' 
shock waves, including p?j? &oc and rj assoc dependences, 
have been investigated in the framework of a hybrid dy- 
namical transport model which consists of two dynamical 
processes, namely parton cascade and hadronic rescater- 
ring. It was found that the splitting parameter D in- 
creases slightly with the transverse momentum of asso- 
ciated hadrons, which results from different interaction- 
lcngths/numbcrs between wave partons and medium in 
strong parton cascade. On the other hand, the pseudo- 
rapidity dependence of splitting parameter D induced by 
a mid-pseudorapidity trigger hadron keeps almost a con- 
stant in mid-pseudorapidity region but drops rapidly in 
higher pseudorapidity region, which is on account of dif- 
ferent strength of jet-medium interactions in the medium 
which has different energy densities in the longitudinal 
direction. Therefore, the correlative researches on the 
properties (j>t and 77 dependences) of 1 Mach-like' corre- 
lation can provide a potential tool to explore the char- 
acters of both partonic and hadronic interactions in the 
formed matter at RHIC. 
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